Dynamic ductile to brittle transition 
in a one-dimensional model of viscoplasticity 
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We study two closely related, nonlinear models of a viscoplastic solid. These models capture essential 
features of plasticity over a wide range of strain rates and applied stresses. They exhibit inelastic 
strain relaxation and steady flow above a well defined yield stress. In this paper, we describe a first 
step in exploring the implications of these models for theories of fracture and related phenomena. 
We consider a one dimensional problem of decohesion from a substrate of a membrane that obeys 
the viscoplastic constitutive equations that we have constructed. We find that, quite generally, 
when the yield stress becomes smaller than some threshold value, the energy required for steady 
decohesion becomes a non-monotonic function of the decohesion speed. As a consequence, steady 
state decohesion at certain speeds becomes unstable. We believe that these results are relevant to 
understanding the ductile to brittle transition as well as fracture stability. 



I. INTRODUCTION 

A wide range of evidence points toward the necessity of 
including plastic deformation near a crack tip among the 
relevant degrees of freedom in theories of dynamic brittle 
fracture. Our own ideas about this issue emerge from 
our recent attempt to study fracture stability using a 
class of cohesive-zone models in which such deformations 
are necessarily absent. As described in our report on 
this project H, we discovered both mathematical and 
physical difficulties that, so far as we could tell, can be 
overcome only by introducing tip blunting or, perhaps, 
a spatially extended process zone in which irreversible 
deformations of the brittle solid are taking place. 

A successful theory of dynamic fracture also must ex- 
lain failure of materials that can flow plastically. Stroh 
understood that cleavage in such materials is an in- 
herently dynamic process in which plastic flow is slow 
enough that stresses can increase to large values near 
the crack tip. In order to explain slower ductile crack 
propagation, McClintock || introduced a novel void nu- 
cleation and coalescence mechanism. At present, theories 
of these two failure mechanisms are separate and incom- 
patible with one another. 

Attempts to explain the differences between ductile to 
brittle fracture usually focus on the emission and mobil- 
ity of dislocations near the crack tip A different 
microscopic mechanism must be responsible for the duc- 
tile to brittle transition in noncrystalline materials such 
as toughened thermoplastics j|j7| . The study of Freund 
and Hutchinson aimed at understanding dynamic frac- 
ture within a theory of continuum plasticity [^|,^). Un- 
der the assumption of elastic strain-rate dominance near 
the tip, they found a non-monotonic dependence of the 
fracture toughness on the crack speed. They used an 
idealized viscoplastic constitutive law, however, in which 
the plastic strain rate is identically zero below a yield 
stress and responds instantaneously to changes in the 
stress. The condition of elastic strain-rate dominance 



breaks down for slow cracks, and thus the analysis of 
Freund and Hutchinson cannot be generally valid. A 
different approach by Freund's [^0) has been to use a 
rate-dependent cohesive-zone model |ll| and to consider 
two separate fracture criteria, one based on stress and 
the other on displacement. This is the class of models 
that we found to be ill-posed for our stability analyses, 
presumably because they omit essential features of an 
extended plastic process zone. 

In a recent study, M. L. Falk and one of the present 
authors (JSL) have proposed a theory of viscoplasticity 
in amorphous materials J12] (denoted FL in what follows) 
in which the equations of linear elastodynamics are sup- 
plemented by nonlinear equations of motion for a set of 
internal state variables. This theory is based directly on 
molecular-dynamics simulations which revealed the ex- 
istence of localized weak regions, called "shear transfor- 
mation zones." The new internal state variables describe 
the populations of these zones, and the nonlinear equa- 
tions describe how those populations determine the time- 
dependent elastic and plastic behavior of the material. 

In the picture presented in FL, the shear transforma- 
tion zones are two-state systems, and transitions between 
those states produce increments in the plastic strain. Be- 
cause a zone that has transformed once cannot transform 
again in the same direction, the plastic strain remains 
bounded when the applied stress is small. An additional 
assumption of FL is that the zones are created and an- 
nihilated at rates determined by the inelastic shear rate. 
It is this process, in which new zones replace old ones, 
that produces persistent plastic flow at sufficiently large 
stresses. Plasticity is a fully dynamic phenomenon in 
this theory. It occurs in a well defined way, depending 
on the state of the system, in response to time dependent 
perturbations. Thus this new theory may be capable of 
describing both brittleness and ductility in fracture. 

We report here on our initial attempts to incorporate 
some of the basic features of the FL theory into a model 
of fracture. As a first step in exploring the implications 
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of this theory, we examine a one-dimensional model of 
decohesion of a thin membrane from a substrate, where 
the membrane obeys a simplified version of the FL vis- 
coplastic constitutive equations. The membrane is pulled 
from the substrate by weak springs, and the point of de- 
cohesion propagates at constant speed, like a crack tip. 
As noted by Barber et. al. [^3) and Marder and Gross 
fL4j , the inverse stiffness of the driving springs in such 
models is analogous, at least in some ways, to the width 
of the strip in two-dimensional mode-I fracture. This 
effective width plays an important role in our interpreta- 
tion of the brittle-ductile transition in this system. The 
one-dimensional nature of the model is, of course, a sig- 
nificant limitation. 

This article is organized as follows. We introduce our 
simplified FL model in Section [o|. In particular, we in- 
troduce two different versions of the nonlinear term that 
determines the plastic yield stress, and we examine the 
time dependence of the system near the onset of plastic 
flow for these two cases. In Section III, we discuss the 
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pi 



linear version of this model and show that it corresponds 
to conventional viscoelasticity with inelastic deformation 
but no persistent plastic flow. Then, in Section IV, we de- 



scribe both analytic and numerical studies of the nonlin- 
ear model, and show how the two rate-dependent mecha- 
nisms introduced in Section [H] produce different versions 
of a brittle-ductile transition. We conclude in Section |v| 
with a discussion of how the lessons learned from this 
simple class of models might be applicable to more real- 
istic situations. 



II. ONE-DIMENSIONAL MODEL OF 
DECOHESION 

We consider a thin membrane decohering from a sub- 
strate as shown in Fig. [jj The system has translational 
symmetry in the direction perpendicular to the plane of 
the Figure. As an additional simplifying approximation, 
we allow the membrane to move only in the direction 
perpendicular to the substrate. The configuration of the 
membrane is thus described by its displacement u(x,t), 
a non-negative function of position x and time t. A co- 
hesive force, 



/ = 







for < u < 5 
otherwise, 



(2.1) 



with a finite range 5, attracts the membrane to the sub- 
strate. Decohesion is driven by weak springs of strength 
a 2 whose relaxed positions are at u = u^. The total 
strain in the membrane is e tot = du/dx. This strain 
and, equivalently, the displacement u consist of additive 
clastic and plastic parts: 



u = u e i + u p i ; 



du e i du p i 



dx dx 
For later notational simplicity, we write 



(2.2) 



dx 



(2.3) 



By definition, the elastic part of the strain is linearly 
proportional to the stress, a — 2p(du e i/dx), where p is 
the shear modulus. Then the equation of motion for the 
membrane is 

d<7 

pu = — k 2 u Q(5 — u) — a 2 (u — Uoo), (2-4) 

ox 

where p is the linear mass density and O(-) is the Heavi- 
side step function. Dots denote time derivatives. 




FIG. 1. A model of one dimensional decohesion driven by 
weak springs. 



The equation of motion for the plastic strain that we 
shall use here is: 
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(2.5) 



This is a simplified version of Eq. (3.14) in FL. (We have 
evaluated the right-hand side of the latter equation in 
the limit of small stress a and have set tia — A, ntot = 
rioo-) By making this small-cr approximation, we lose 
some of the memory effects that were obtained in FL 
via a strongly nonlinear c-dependence of the rate factors 
in the equation for e. We believe that the absence of 
those effects makes only quantitative and not qualitative 
differences in the results to be presented here. However, 
that point may require further investigation. 

Our single state variable, A(x,t), is a measure of the 
imbalance in the populations of the two- stat e systems. 
The first term on the right-hand side of ([2~5|), A<t/t, is 
the usual linear relation between the plastic strain rate 
and the stress. The second term, — A/r, is the rate at 
which these two-state systems transform spontaneously 
from their "forward" to their "backward" states, and is 
therefore a negative contribution to e. 

Our equation of motion for A has the form: 



A = e- F(e,a)A. 



(2.6) 



This is exactly the same as Eq. (3.15) in FL except that 
we have not yet specified the strain-rate dependent cou- 
pling T between A and A. The first term on the right 
hand side of (2.6) simply expresses our assumption that 
the transitions within the two-state systems correspond 
to increments in plastic strain. The second, i.e. the non- 
linear term, is the effect of creation and annihilation of 
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these zones and, as we shall see immediately, is respon- 
sible for the existence of a finite plastic yield stress. 

To illustrate the properties of the nonlinear term in 
(2.6), we consider two plausible forms of T that pro- 
duce qualitatively different dynamic decohesion in certain 
regimes. The form of T is restricted by the assumption 
that it must vanish when the plastic strain rate vanishes, 
ft is thus proportional to some power of e. (In higher 
dimensions, we would also require rotational invariance.) 
The first form that we shall examine is the same as that 
used in FL: T\ = jiea. Here, the coupling T\ is pro- 
portional to the local rate of plastic energy dissipation. 
We call this Model 1. Note that T\ can be negative in 
some circumstances. Our second possibility, Model 2, is 
one in which only the local plastic strain rate controls 
the evolution of A, in which case the simplest choice is 

?2 = ll<?- 




FIG. 2. The steady state values of A in Model 1 as func- 
tions of ct. The two curves Aia = Act and Aib = I/71C cross 
at (ijyi, A y i). The arrows indicate the direction of motion of 
A for fixed ct. 



This rate vanishes when the stress approaches yield from 
above. However, the relaxation time n for perturbations 
away from the flowing state diverges for ct near a v \ : 



n 



(2.11) 



'yil 



This is quite unlike the conventional clastic-idcally plastic 
solid in which this relaxation time is zero by definition. 

A possibly unphysical feature of this model is that, for 
some initial conditions, A may increase i ndefi nitely as a 
function of time. It is easy to see from ( |2.7| ), however, 
that if A < A y i at any time, it will remain so for all 
other times regardless of the stress history. 




FIG. 3. Steady state values of A as functions of ct for Model 
2. One of the three, A2B+, is unstable to small perturbations. 
Since the steady state curves never cross, the system remains 
close to the non-flowing steady state A2A if the stress in- 
creases slowly enough. The meaning of the arrows is the same 
as in Fig. 



We explore first the beh avio r of M ode l 1. Substituting 
the expression for e from (|2.5|) into (2.6), we obtain 
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(Act - A) (1 - 7 icrA) . 



(2.7) 



For constant stress ct, there are two stationary solutions: 

1 



A = Aia = Act, and A = A 1B 
which, as shown in Figure ||, cross at 

A = Ayi = 



71 CT 



1 



CT = CT„i 



/71A 



(2.8) 



(2.9) 



At any fixed ct, only one of these stationary solutions is 
stable against perturbations. For a < a y i, the stationary 
solution Aia with e = is stable. For ct > a y i, on the 
other hand, the stable stationary solution A = Axb is a 
flowing steady state with 



CTT a 



(2.10) 



Model 2 exhibits an important qualitative difference in 
its behavior. Let us perform the analy sis of the p rece ding 
paragraphs using T 2 . Substituting ( pH| ) into (2.6), we 
obtain 



A = — (Act — A) (^a 2 -^A + 1 

T \ T T 



(2.12) 



We again look for the stationary states A = 0. The 
situation is shown in Fig. [| In this case, the state with 
A = A 2 a = Act and no plastic flow, e = 0, is stable for all 
ct. As seen in Fig. [j[ it never intersects any other state 
in the A-ct plane. For 
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A V 72 

two new stationary states appear with 
A 



A = A 2B ± = - [o ± \J<J 2 - <J y2 



(2.13) 



(2.14) 



The state with A = A2B- is the stable one of the pair. 
The plastic strain rate in this stationary state is non-zero: 



3 



>l,2 



(steady state, a > <J y -z). 

(2.15) 



The characteristic decay time ti of perturbations around 
the flowing stationary state again diverges at the yield 
stress, although the divergence is not as strong as in 
Model 1: 
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(a > a 



V2) 



'y2 



(2.16) 



The two nonlinear models exhibit a number of sim- 
ilar features. Most importantly, steady plastic flow in 
response to a stress above a yield level is a natural con- 
sequence of the dynamical constitutive equations. The 
flow has a non-zero response time to changes in the stress. 
It can also be shown that inelastic strain is partially re- 
covered in both models. However, there are several im- 
portant differences between the two nonlinear models. 
First, the steady flow rate does not vanish in Model 2 at 
c = Cy2- Second, "runaway" behavior cannot occur in 
that model since A < for A > A2A- And finally, for 
stresses greater than the yield stress, there are two stable 
stationary solutions in Model 2 as opposed to only one 
in Model 1. Which one of these is selected by the system 
depends on the stress history. For example, only the non- 
flowing state in Model 2, with A = A 2 a, occurs if the 
stress is increased slowly enough, i.e. when &/a <C 1/r. 
As we shall see, this distinction between the models leads 
to qualitatively different behaviors at small speeds. 

It is convenient at this point to convert to dimension- 
less variables in which all lengths are measured in units 
of the range of the cohesive interaction <5, time is mea- 
sured in units of 5^J p/2/z, and stress in units of 2/i. For 
simplicity, we continue to use the symbols u and a for 
our dimensionless displacements and stresses. We also 
restrict our attention to steady-state solutions moving in 
the negative x direction with speed v. All functions of x 
and t in the frame of reference moving with the decohe- 
sion front depend only on the combination x' = x + vt. 
Without loss of generality, we set x' = at the point 
of decohesion where, in these units, u = 1. Then, for 
simplicity, we set x' = x. 

Our equations of motion now have the form: 



v 2 (u"i + Up/) = u'et ~ « 2 K* + u p i)0(-x)- 



' P i 
vA' 



a (u ei + u p i - Uoo), 
- (Hi ~ A) , 

T 

vu'pi -^{vu'pi,u' el )A, 



(2.17) 
(2.18) 
(2.19) 



where primes denote differentiation with respect to x, 
and 



T 2 = l2(vupi) 2 . 



(2.20) 
(2.21) 



Finally, we derive an expression for the decohesion re- 
sistance G(v), which is the work that the driving springs 
perform on the membrane per unit length of advance of 
the decohesion front. Since the driving springs relax to 
their equilibrium length far behind the decohesion front, 
all of their stored elastic energy ahead of the front must 
be dissipated in the decohesion process. Thus, the total 
work done must be 
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G(v) — —a u 



2„,2 

00 ■ 



(2.22) 



If we multiply ( |2.17| ) by u'(x) and integrate over x, we 
obtain 



G{v)= l -^ + 



dx 



l pl u 'el + 2 ( U f>') 2 



(2.23) 



The first term on the right-hand side of ( 2.23 ) is clearly 
the energy spent in breaking cohesive bonds. The sec- 
ond term is the energy dissipated in plastic work in the 
neighborhood of the decohesion front. The third term is 
the energy locked into the plastic wake left by the deco- 
hesion. The decohesion front in this model may leave a 
residual plastic strain behind it. Our problem is to com- 



pute G{v) explicitly as a function of v and then use (2.22) 
to determine v as a function of the driving force auoo . 



III. LINEAR ANALYSIS 



Before going ahead with an analysis of this nonlinear 
model of viscoplasticity, it will be useful to look briefly 

at the linear case in which we set T 

„,/ 



we have A 
is 



in (2.6). Then 



pi 



; and the equation of motion for u p i 



Upl 



pi 



= -(Au e 



hi) 



(3.1) 



The remaining equation of motion is (2.4) or, equiva- 
lently, @). 

Far away from the region where decohesion is taking 
place, our system is translationally invariant, and we can 
compute a dispersion relation for waves of the form u = 
uoexp(ikx — iuJkt). In the limit of vanishing a, we find 



The wave speed c is 



A 



1 - lUJ k T 



1 



c = lim Re — ^ = — = . 

k^o k Vl + A 



(3.2) 



(3.3) 



It is important to recognize that, by linearizing, we 
have reduced our system to a conventional model of vis- 
coelasticity. The solution of the time-dependent version 
of (3.1) can be written in the familiar form 



etot(t) = ^a(t)-X 



/dt' exp 
-00 



1 



(t-0 



«r(f) (3.4) 
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where e iot is the total (elastic plus plastic) strain, and a = 
du e i/dx is the stress in dimensionless units. Equivalently, 



a{t) =c 2 6 to t{t) + \ \ dt'exp 



1 



-Z3(*-0 



TCT 



e t ot(t'). 
(3.5) 



From (3.4) we find that the creep compliance — the vari- 
ation of the strain that is produced by a unit jump in the 
stress — is 



C(i) = l + A(l-e- t/r 



(3.6) 



The system exhibits unit instantaneous elasticity, follow- 
ing which the strain increases on the time scal e r to its 
final value 1 + A = 1/c 2 . Similarly, we see from (3.5) that 
a unit jump in the strain produces first an instantaneous 
jump in the stress, after which the stress decreases to a 
constant, nonzero value. 



Because our equations of motion (2.17) and (3.1) are 
linear, we can solve the decohesion problem analytically. 
(We shall need these linear solutions in order to interpret 
features of the nonlinear solutions described in the next 
Section.) The analysis is particularly simple if we take 
the limit of weak driving springs, a — > 0. To do this, 
we must also take the limit — > oo in such a way as 
to keep auoo constant. That is, we keep G(v) fixed in 



(2.22) 



For x < 0, we can set a = immediately in (2.17) and 
look for a solution in the form: 



u e i = Ae qx , 



u p i 



Be qx . 



(3.7) 



That is, we look for values of A, B, and q that satisfy 
the homogeneous equation: 



-/3V 

-X/VT, 



k 2 , v 2 q 2 + k 2 



1/VT 



0. 



The solvability condition for (3 

q 3 VTj3 2 +q 2 f3 2 - qvrn 



is: 



.2 



0. 



(3.8) 



(3.9) 



(la 



a 



12 = 



+ 0(a 2 ). 



vt/3 2 

We therefore construct solutions of the form: 



u el = D ie qoX + D 2 e^ x + c 2 Uoo , 



XD 1 



1 + vrq 

(1 - C^Uoc 



„qox 



XDo 



1 + vrq 2 



(3.12) 



(3.13) 



(3.14) 



Here, we have included the particular solutions (simple 
constants) that satisfy the boundary condition u — > Uoc 
at x — * +oo. We then require that u e i and u p i and their 
derivatives be continuous at x = 0. Calculating to first 
order in aiioa, we find: 



D\ = —c u 
Do 



A^+c 2 u c 



'12 



9o 

% 

) 

92 



A 1- 



<12 



where A is given in (3.11). Finally, 



K(v) = y/2G(v) = au x = 



qic ( (3 2 VTq 1 + (3 2 



0c \ c 2 vrqi + 1 



(3.15) 
(3.16) 

(3.17) 



K(v) 



0.5 1 

v/c 

FIG. 4. Decohesion toughness K(v) for the linear model 
with t = 10 and A = 1. 



where 



B 2 



1 



01 EE 1 - V 2 /C 2 



(3.10) 



Eq. ( |3.9D has only one positive root, say q±. The bound- 
ary condition u e i (0) + u p i (0) = 1 is therefore sufficient to 
determine uniquely the solution in the region x < 0. We 
find: 



1 + qiVT 



1 



qivrc* 



B 



Ac 2 



1 + qivrc 



(3.11) 



For x > 0, all k's appearing in ( |3.9| ) must be replaced 
by a's. The resulting equation has two negative solutions 
for q which, for small a, are: 



We show a representative graph of if as a function 
of the speed v in Fig. ||. As expected from a model of 
this kind, K(v) is monotonic, diverges at v = c as 1//3 C , 
and is equal to k for v = 0, which confirms that viscous 
dissipation is negligible for slow decohesion. 



IV. SOLUTIONS OF THE NONLINEAR MODELS 

We turn now to the nonlinear models defined by 
Eqs. (2.17) through ( [2.19] ). These are equivalent to a 
set of ordinary differential equations which, for Model 1, 
are: 
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/?V = - (Act - A) + k 2 @(-x) u + a 2 (u - u M ), (4.1) 

T 

e' = — (A<7-A), (4.2) 

A' = — (Xa - A) [1 - 7 iaA] , (4.3) 



it 



cr + e. 



(4.4) 



To obtain the analogous equat ions for Model 2, replace 
7i<r in the square brackets in (4.3) by 72«e' = 72(Acr — 
A)/r. ^ 

We integrate these equations numerically using a 
twelfth-order predictor-corrector algorithm. The initial 
conditions are CT = e = A = 0atx— * — oo, far ahead of 
the decohesion front. Our strategy is to fix the material 
parameters A, 7, and r (or, equivalently, c, <r y , and t), 
the velocity v, and the strength of the driving springs a 2 , 
and to adjust ttoo to obtain a solution with the property 
that u(0) = 1 and u — > Uqo as x — > +00. Such a solu- 
tion always exists. Recall that, in our analogy with the 
crack propagating in a prestressed strip, the parameter 
Moo is analogous to the displacement of the strip edges 
far ahead of the crack tip. We therefore might think of 
our procedure as adjusting the driving stress on the strip 
to achieve a certain velocity of fracture propagation. 

Before looking in detail at these solutions, consider 
the following thought experiment. Imagine that we start 
with a static, unstressed system and — 0. Suppose 
also that the cohesive springs act only for a; < 0; that 
is, we arbitrarily disconnect them for x > 0. Let us now 
increase Uoo from zero quasistatically. In this limit of in- 
finitesimally slow displacement and a fixed position of the 
decohesion front, the nonlinear models are indistinguish- 
able from the linear model as long as the stress in the 
membrane nowhere exceeds the plastic yield stress. This 
is because, for a < a y , the quasistatic system must stay 
arbitrarily close to the nonflowing state with e = A = Xa 
and the nonlinear term in (^^) is irrelevant. 

The linear theory tells us that the largest stress occurs 
at x — where, for this quasistatic situation, it has the 
value Umax = ctUooK. At this point, the displacement of 
the membrane is u(0) — otUoo/c. Clearly, the behavior of 
this system depends sensitively on whether or not a max 
exceeds the plastic yield stress a y before u(0) reaches the 
breaking point u(0) = 1. Thus the critical value of the 
yield stress that marks some kind of quasistatic boundary 
between brittle and ductile behavior of these models is 

(Ty = KC. 

If <7y > kc, then the threshold for propagating deco- 
hesion is reached before any plastic flow occurs, and we 
deduce that both nonlinear models behave much like the 
linear model for small enough speeds v — i.e. they are 
"brittle." On the other hand, if a y < kc, then plastic flow 
occurs before the leading cohesive spring breaks. In this 
case, the two nonlinear models behave differently from 
one another. 

In Model 1, plastic flow must begin as soon as the 
maximum stress reaches the yield stress. As shown in 



Fig. ||, the flowing and nonflowing states cross at this 
point. If is increased very slowly beyond this point, 
just as in conventional models of plasticity, the stress at 
x = remains fixed at a y . The displacement u(0) also 
remains fixed at its value below the decohesion threshold, 
u(0) — 1, because no additional stress can be applied to 
stretch the cohesive springs. The only thing that can 
happen is that, as continues to grow, the material in 
the region x > deforms plastically. Thus, decohesion is 
not initiated, but an indefinitely large amount of plastic 
work is done on the system. We therefore anticipate that 
the decohesion toughness K(v) for Model 1 must diverge 
at v — whenever a y < kc. 

The most interesting questions, of course, have to do 
with the behavior at nonzero propagation speeds v, where 
the quasistatic approximations are not necessarily valid. 
In general, our decohesion criterion u(0) = 1 implies that 
the breaking stress <j(0) increases with increasing v. (We 
shall not consider a stress-based criterion, which might be 
simpler in some respects.) Actually, at nonzero speeds, 
the stress reaches its maximum some distance behind the 
decohesion front. In the linear version of the model, this 
maximum stress diverges at v = c. Thus, even if the 
maximum stress kc is less than a y at v — 0, it will become 
greater than a y at some onset speed for plastic flow that 
we shall call v p . At speeds larger than v p , the system 
must deform plastically. 




FIG. 5. Decohesion toughness K(v) with a c — 3kc and 
r = 10, for three values of the driving a. The limit of weak 
driving a — » exists only for v < v p w 0.73c. 



To see what happens at v p , we m ust look at th e nu- 
merical solutions of Equations (4.1) through (4.3). We 
continue to consider only Model 1 for the moment. In 
Fig. ||, we show the decohesion toughness K{v) — au c 



as a function of v/c, for a y = 3kc, r = 10, A = 1, and for 
three different values of a. For comparison, we also show 
K(v) for the linear theory. The most striking feature 
is that these four curves are almost coincident for v/c 
less than a critical value of about 0.73; but they break 
away from the linear theory at larger speeds, the systems 
with smaller a being the most dissipative. We see even 
more explicitly what is happening in Figs. (^|) and (|7]) 



G 



where we have plotted the total displacement u(x) and 
the plastic strain e(x) for a = 10~ 3 , and for two veloci- 
ties v/c = 0.71 and 0.75, just below and just above the 
critical speed respectively. Note that the plastic strain is 
very much larger for the slightly larger value of v. In that 
case, the plastic strain grows almost linearly with x be- 
fore it reaches its peak. The spatial extent of the region 
in which this plastic strain accumulation happens seems 
to scale linearly with a -1 . Note also that decohesion 
produces a residual plastic deformation and, accordingly, 
a residual stress (not shown here) that persist infinitely 
far behind the front. 



all (x) 



v/c = 0.71 
v/c — 0.75 




FIG. 6. Total displacement U(x) as a function of x, both 
measured in units of a" 1 , for speeds just below and just above 
v p . A considerable increase in the driving force alloc is needed 
to increase the decohesion speed by 4%. 
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FIG. 7. Plastic strain e(x) as a function of ax for speeds 
just above and just below v p . Note that the peak in e(x) for 



v/c 



0.75 > 



is more than twice that of the peak for 



v/c = 0.71 < v p , and the level of plastic strain far behind the 
decohesion front changes by a factor of about 5. 



v p is the speed at which the maximum stress just exceeds 
the plastic yield stress. To confirm this interpretation, in 
Fig. ^ we plot v p as a function of a y and compare this 
with the prediction of the linear model for the maximum 

stress 



<jy at v — Vp(ay), (4-5) 



where aitoo, is given in Eq. ( 3.17 ). Agreement with the 
linear theory becomes exact in the limit of slow decohe- 
sion. We can also qualitatively understand the fact that 
the prediction of the linear model consistently overesti- 
mates the onset velocity v p . When the nonlinear term 
can be treated perturbatively, it leads to a decrease in 
the relaxation rate of A, since the right hand side of 
Eq. ( p~3| ) is reduced. The right hand side of the equation 
for the derivative of the stress (4.1) is therefore increased 
thus allowing the stress to reach a higher level before the 



a (u 



3 ) term reverses the sign of a 1 . As a result, the 



maximum stress in a nonlinear system reaches yield at a 
lower velocity than in the corresponding linear system. 



o.o 



0.4 




Model 1 O 
Prediction of the linear model 
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a yl /KC 

FIG. 8. Onset velocity v v as a function of the yield stress 
for Model 1. 



So far, we have looked in detail only at the physically 
less realistic situation in which the plastic yield stress 
is higher than the breaking stress at the v — threshold 
for decohesion. We now consider the case where the yield 
stress is lower than the breaking stress. In Model 1, if 
cty < kc, then we are always in the regime where the 
a — > limit fails to exist. In examining the behavior in 
this regime, therefore, we choose a small, fixed value for 
a, specifically a = 0.01, and look at various values of 
other parameters. 



We deduce from this data that the a — > limit exists 
only for v < v p . That is, if we try to drive decohesion 
at a speed greater than v p with springs of arbitrarily 
small force-constant a 2 , the dissipated energy per unit 
advance of the front grows without bound. For Model 1, 
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0.5 1 

v/c 

FIG. 9. Decohesion toughness K(v) for Model 1 in the case 
a y i = 0.5kc for q = 0.01. 



Fig. H is a graph of K(v) as a function of v/c for 
(jy = 0.5kc. The most interesting feature of this graph 
is that dK/dv < for speeds v between zero and, say, 
v a . Propagation at speeds in that interval must be un- 
stable; thus, if we increase the driving parameter to 
some value such that K = auoo > K(v a ), then v must 
jump to some value on the high-speed, stable branch of 
this curve. Conversely, if we decrease the driving force so 
that craoo < K(v a ), then decohesion must stop abruptly. 
As anticipated, the decohesion toughness is large at small 
speeds because the plastic strain relaxes very slowly near 
threshold, and the flowing region extends far behind the 
decohesion front. At larger speeds, the deformation is 
more localized, and there is less dissipation. At yet larger 
speeds, of course, the driving force must increase in or- 
der to make the front move at speeds comparable to c. 
To illustrate these variations in the plastic deformation 
explicitly, in Fig. |l0| we plot the plastic strain rate e as 
a function of x/v. This figure can be interpreted as the 
plastic strain rate as a function of time after passage of 
the decohesion front. 
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FIG. 10. Plastic strain rate as function of time x/v after the 
passage of the decohesion front in Model 1. The parameters 
are the same as in Figure |^. Plastic flow persists much longer 
for slow decohesion. Also note that plastic strain recovery is 
appreciable only when decohesion is fast. 



We turn now to the properties of Model 2, which must 
behave in a less conventional manner according to Fig. |^. 
Even when the breaking stress exceeds the plastic yield 
stress, a y < kc, the system can remain in the station- 
ary state with A = A2A, e = so long as the stress is 
raised sufficiently slowly. No plastic flow occurs, and the 
decohesion toughness at v — > must be k, just as in the 
linear model. 

To predict the onset of plastic flow in Model 2, i.e. to 
compute the analog of v p (a y ), we can use the linear the- 
ory to estimate when the nonlinear term in the equation 
of motion ( |4.3D becomes non-negligible. It is useful to 
carry out this exercise for both models. For Model 1, 
validity of the linear approximation requires 



71 erA = liu' el u' pl < 1- 
Our linear analysis tells us that 



7iA(acUoo)' 



1 - 



vrqi 
1 

1 + vrqi 



-x/vT 



(4.6) 



(4.7) 



This quantity vanishes at the decohesion front, x = 0, 
and rises monotonically to a constant as x — > 00. We 
know that 71A = 1/cr 2 , and, for kvt <C 1, aitm ~ n. Thus 
the inequality in ([O]) reduces to a y i kc, consistent 
with our quasistatic analysis for this model. 

Mode l 2 behaves differently. The analog of the inequal- 
ity (111) is 



7 2 w'A S l2vu' pl u'^ < 1. 
From the linear analysis, we find 

(1 + vrqi) 



(4.8) 



l2vu' pl u' pl w j2vX 2 qfc 4 



1 



(1 
vrqi 



vrqic 2 ) 2 

„~X/VT 



e- x/vT x 



1 



vrqi 



(4.9) 



Now the nonlinear correction is localized in a finite region 
whose size is of order vt near the decohesion front. We 
use 72 A 2 = 4t/cj 2 , 9 . Then, in the case kvt <C 1, the 



inequality (4.8) becomes 



vt <C 



or, equivalcntly, 



KVT 



u y2 
Ack 3 ' 



<Jy 2 



2a(0) 



(4.10) 



(4.11) 



Both sides of ( 4.11 ) are ac cur ate only for kvt <C 1. In the 
opposite limit, kvt 3> 1, (4.8) reduces simply to a V 2 k. 
The important point is that, when v is sufficiently small 
in Model 2, brittle behavior can occur for values of the 
plastic yield stress a y 2 that are smaller than the decohe- 



sion stress <r(0). The right-hand side of (4.10) gives us 
an upper bound for v p t for Model 2. 



8 



15 



10 



K(v) 







10~ 3 — 






1Q -4 ___ 




a 2 = 


10 -5 .... 

















v/c 



().: 



FIG. 11. Decohesion toughness K(v) in Model 2 for three 
different values of a. Other parameters are a y 2 = 0.8kc, A = 1 
and t — 4. 



These features of the behavior of Model 2 are confirmed 
by our numerical results. In Fig. [ll], we show the deco- 
hesion toughness K(v) as a function of v for a y 2 — 0.8nc 
for three different values of a 2 , A = 1 and r = 4. As 
expected, K(0) = k. There is a stable region at small v 
where dK/dv > and in which the limit a — > exists. 
The the onset of plastic flow speed, and the failure of 
that limit, occurs at v = v p where K(y) first begins to 
rise sharply. Beyond its maximum, the fracture tough- 
ness in Model 2 behaves qualitatively like Model 1. That 
is, for small but nonzero values of a, there is an unstable 
region in which K(v) decreases for increasing v. At yet 
larger values of v, K(v) rises again and diverges as v ap- 
proaches c in the limit of a — > 0. As shown in Fig. fl2L the 
plastic flow onset velocity v p in Model 2 vanishes only as 
the yield stress <j v is reduced to zero. Its behavior for 
small yield stresses is consis tent w ith the prediction of 
the perturbation theory Eq. ( 4.10 ). 
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FIG. 12. Onset velocity v p in Model 2 for A = r = 1 as 
a function of the yield stress a y 2- Lower yield stresses are 
numerically difficult to access. 



V. DISCUSSION 

We have explored these one-dimensional models of de- 
cohesion primarily as an attempt to understand how a 
fully dynamic description of viscoplasticity might play a 
role in theories of dynamic fracture. We are particularly 
interested in how the concepts of brittleness and ductility 
will emerge in such theories. 

To begin this discussion, consider one way in which we 
might expect the brittle-ductile transition to appear in 
a theory of, say, mode I fracture along the centerline of 
an infinitely long, two-dimensional strip. Let the width 
of the strip be 2W, and suppose that the driving force is 
produced by a rigid displacement of size Woo at its edges. 
In the case of brittle fracture, the stress-intensity factor 
at the tip of the crack is proportional to Kj = Uoo/y/W. 
Infinitely far behind the crack tip, the steady-state crack 
opening displacement is 21*00. In the limit W — + 00, Kj 
remains fixed for a fixed speed of crack propagation, and 
thus the ratio u^jW vanishes. The crack remains sharp 
and narrow on the macroscopic scale W . 

The extreme ductile version of this situation is one in 
which the solid is replaced by a viscous fluid, and the 
crack becomes a finger in a Hele-Shaw cell. In this case, 
the steady-state finger has a width of approximately 
W, and the dissipation rate, proportional to Kj , diverges 
as W — » 00. In any real two or three dimensional solid, of 
course, the plastic yield stress is nonzero. Therefore, as 
we increase W at fixed crack speed — no matter whether 
the crack is advancing in a brittle or ductile manner — 
we must eventually get to the point where the stresses 
far away from the crack tip drop below the plastic yield 
stress. The dissipation may become very large, but it 
remains finite when W — > 00. 

In our one-dimensional model, the closest available 
analog of the width W is the length scale a -1 . How- 
ever, we have no analog of the stress concentration that 
is produced by a real second dimension, and thus we have 
no way in which the far-field stresses can be made arbi- 
trarily small — less than the plastic yield stress — by 
taking the limit of a large system. In its brittle mode, as 
we have seen, the analog of the stress-intensity factor for 
our decohering membrane is K(y) = au^ ss Uqo/W; this 
quantity remains finite at fixed v in the limit a — > 0. In 
its ductile mode, however, our system is behaving more 
like a finger in a Hele-Shaw cell than a crack in a solid 
strip. As soon as plastic flow starts at any point, the dis- 
sipation rate diverges as a — > 0. In short, the distinction 
between brittle and ductile failure in this system must 
be qualitatively unlike that which occurs in real fracture. 
To understand the latter, we shall have to carry out fully 
two-dimensional investigations. 

What, then, are the lessons to be learned from this 
exercise? What questions does it raise? We have con- 
firmed, as expected, that the FL model of viscoplastic- 
ity produces both brittle and ductile propagating failure 
modes. In this one-dimensional version, the transition 
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between brittle and ductile behavior is perfectly sharp 
and well defined; it is distinguished by the divergence 
of the decohesion toughness in the limit a — > 0. Is there 
any such sharp distinction in higher dimensions? The be- 
havior of our one-dimensional models, especially at small 
propagation speeds v, is strongly sensitive to our choice 
of the nonlinear coupling between the plastic strain rate 
and the new variable A that describes the internal state 
of the system. It will be important to learn whether more 
realistic models in higher dimensions exhibit similar sen- 
sitivity to details of the mechanisms that control plastic 
flow. 

Perhaps the most interesting but problematic aspect 
of our results is related to that sensitivity. Both of the 
models discussed here exhibit unstable steady-state so- 
lutions at low propagation speeds. We know that these 
solutions are unstable because the decohesion toughness 
K (v) decreases with increasing v. In both models, K (v) 
rises again at higher speeds, and the stable high-w so- 
lutions are ductile. In Model 2, however, there is also 
a stable small-t; solution that is brittle. That is, there 
exists a range of values of K(v) for which there are two 
stable steady-state solutions, one brittle and one duc- 
tile. Within such a range of driving forces, the system 
is likely to exhibit complex, non-steady behavior. What 
might be the analog of such behavior in two-dimensional 
models of fracture? Might there be situations in which 
multiple solutions exist but both are brittle? Or might 
the slow solution be the ductile one? More generally, 
might the new dynamics emerging from the FL model of 
viscoplasticity be a clue for understanding the complex 
instabilities and different modes of behavior observed in 
real fracture? 
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